Epigenetic regulators controlling osteogenic lineage commitment and bone formation.

Bone formation and homeostasis are controlled by environmental factors and endocrine regulatory cues that initiate intracellular signaling pathways capable of modulating gene expression in the nucleus. Bone-related gene expression is controlled by nucleosome-based chromatin architecture that limits the accessibility of lineage-specific gene regulatory DNA sequences and sequence-specific transcription factors. From a developmental perspective, bone-specific gene expression must be suppressed during the early stages of embryogenesis to prevent the premature mineralization of skeletal elements during fetal growth in utero. Hence, bone formation is initially inhibited by gene suppressive epigenetic regulators, while other epigenetic regulators actively support osteoblast differentiation. Prominent epigenetic regulators that stimulate or attenuate osteogenesis include lysine methyl transferases (e.g., EZH2, SMYD2, SUV420H2), lysine deacetylases (e.g., HDAC1, HDAC3, HDAC4, HDAC7, SIRT1, SIRT3), arginine methyl transferases (e.g., PRMT1, PRMT4/CARM1, PRMT5), dioxygenases (e.g., TET2), bromodomain proteins (e.g., BRD2, BRD4) and chromodomain proteins (e.g., CBX1, CBX2, CBX5). This narrative review provides a broad overview of the covalent modifications of DNA and histone proteins that involve hundreds of enzymes that add, read, or delete these epigenetic modifications that are relevant for self-renewal and differentiation of mesenchymal stem cells, skeletal stem cells and osteoblasts during osteogenesis.

LINC01638 sustains human mesenchymal stem cell self-renewal and competency for osteogenic cell fate.

The skeleton forms from multipotent human mesenchymal stem cells (hMSCs) competent to commit to specific lineages. Long noncoding RNAs (lncRNAs) have been identified as key epigenetic regulators of tissue development. However, regulation of osteogenesis by lncRNAs as mediators of commitment to the bone phenotype is largely unexplored. We focused on LINC01638, which is highly expressed in hMSCs and has been studied in cancers, but not in regulating osteogenesis. We demonstrated that LINC01638 promotes initiation of the osteoblast phenotype. Our findings reveal that LINC01638 is present at low levels during the induction of osteoblast differentiation. CRISPRi knockdown of LINC01638 in MSCs prevents osteogenesis and alkaline phosphatase expression, inhibiting osteoblast differentiation. This resulted in decreased MSC growth rate, accompanied by double-strand breaks, DNA damage, and cell senescence. Transcriptome profiling of control and LINC01638-depleted hMSCs identified > 2000 differentially expressed mRNAs related to cell cycle, cell division, spindle formation, DNA repair, and osteogenesis. Using ChIRP-qPCR, molecular mechanisms of chromatin interactions revealed the LINC01638 locus (Chr 22) includes many lncRNAs and bone-related genes. These novel findings identify the obligatory role for LINC01638 to sustain MSC pluripotency regulating osteoblast commitment and growth, as well as for physiological remodeling of bone tissue.

Protein arginine methyltransferases PRMT1, PRMT4/CARM1 and PRMT5 have distinct functions in control of osteoblast differentiation.

Osteogenic differentiation of mesenchymal cells is controlled by epigenetic enzymes that regulate post-translational modifications of histones. Compared to acetyl or methyltransferases, the physiological functions of protein arginine methyltransferases (PRMTs) in osteoblast differentiation remain minimally understood. Therefore, we surveyed the expression and function of all nine mammalian PRMT members during osteoblast differentiation. RNA-seq gene expression profiling shows that Prmt1, Prmt4/Carm1 and Prmt5 represent the most prominently expressed PRMT subtypes in mouse calvarial bone and MC3T3 osteoblasts as well as human musculoskeletal tissues and mesenchymal stromal cells (MSCs). Based on effects of siRNA depletion, it appears that PRMT members have different functional effects: (i) loss of Prmt1 stimulates and (ii) loss of Prmt5 decreases calcium deposition of mouse MC3T3 osteoblasts, while (iii) loss of Carm1 is inconsequential for calcium deposition. Decreased Prmt5 suppresses expression of multiple genes involved in mineralization (e.g., Alpl, Ibsp, Phospho1) consistent with a positive role in osteogenesis. Depletion of Prmt1, Carm1 and Prmt5 has intricate but modest time-dependent effects on the expression of a panel of osteoblast differentiation and proliferation markers but does not change mRNA levels for select epigenetic regulators (e.g., Ezh1, Ezh2, Brd2 and Brd4). Treatment with the Class I PRMT inhibitor GSK715 enhances extracellular matrix mineralization of MC3T3 cells, while blocking formation of H3R17me2a but not H4R3me2a marks. In sum, Prmt1, Carm1 and Prmt5 have distinct biological roles during osteoblast differentiation, and different types histone H3 and H4 arginine methylation may contribute to the chromatin landscape during osteoblast differentiation.

Identification of molecularly unique tumor-associated mesenchymal stromal cells in breast cancer patients.

The tumor microenvironment is a complex mixture of cell types that bi-directionally interact and influence tumor initiation, progression, recurrence, and patient survival. Mesenchymal stromal cells (MSCs) of the tumor microenvironment engage in crosstalk with cancer cells to mediate epigenetic control of gene expression. We identified CD90+ MSCs residing in the tumor microenvironment of patients with invasive breast cancer that exhibit a unique gene expression signature. Single-cell transcriptional analysis of these MSCs in tumor-associated stroma identified a distinct subpopulation characterized by increased expression of genes functionally related to extracellular matrix signaling. Blocking the TGFß pathway reveals that these cells directly contribute to cancer cell proliferation. Our findings provide novel insight into communication between breast cancer cells and MSCs that are consistent with an epithelial to mesenchymal transition and acquisition of competency for compromised control of proliferation, mobility, motility, and phenotype.

Epigenetic and transcriptome responsiveness to ER modulation by tissue selective estrogen complexes in breast epithelial and breast cancer cells.

Selective estrogen receptor modulators (SERMs), including the SERM/SERD bazedoxifene (BZA), are used to treat postmenopausal osteoporosis and may reduce breast cancer (BCa) risk. One of the most persistent unresolved questions regarding menopausal hormone therapy is compromised control of proliferation and phenotype because of short- or long-term administration of mixed-function estrogen receptor (ER) ligands. To gain insight into epigenetic effectors of the transcriptomes of hormone and BZA-treated BCa cells, we evaluated a panel of histone modifications. The impact of short-term hormone treatment and BZA on gene expression and genome-wide epigenetic profiles was examined in ERαneg mammary epithelial cells (MCF10A) and ERα+ luminal breast cancer cells (MCF7). We tested individual components and combinations of 17ß-estradiol (E2), estrogen compounds (EC10) and BZA. RNA-seq for gene expression and ChIP-seq for active (H3K4me3, H3K4ac, H3K27ac) and repressive (H3K27me3) histone modifications were performed. Our results show that the combination of BZA with E2 or EC10 reduces estrogen-mediated patterns of histone modifications and gene expression in MCF-7ERα+ cells. In contrast, BZA has minimal effects on these parameters in MCF10A mammary epithelial cells. BZA-induced changes in histone modifications in MCF7 cells are characterized by altered H3K4ac patterns, with changes at distal enhancers of ERα-target genes and at promoters of non-ERα bound proliferation-related genes. Notably, the ERα target gene GREB1 is the most sensitive to BZA treatment. Our findings provide direct mechanistic-based evidence that BZA induces epigenetic changes in E2 and EC10 mediated control of ERα regulatory programs to target distinctive proliferation gene pathways that restrain the potential for breast cancer development.

LncMIR181A1HG is a novel chromatin-bound epigenetic suppressor of early stage osteogenic lineage commitment.

Bone formation requires osteogenic differentiation of multipotent mesenchymal stromal cells (MSCs) and lineage progression of committed osteoblast precursors. Osteogenic phenotype commitment is epigenetically controlled by genomic (chromatin) and non-genomic (non-coding RNA) mechanisms. Control of osteogenesis by long non-coding RNAs remains a largely unexplored molecular frontier. Here, we performed comprehensive transcriptome analysis at early stages of osteogenic cell fate determination in human MSCs, focusing on expression of lncRNAs. We identified a chromatin-bound lncRNA (MIR181A1HG) that is highly expressed in self-renewing MSCs. MIR181A1HG is down-regulated when MSCs become osteogenic lineage committed and is retained during adipogenic differentiation, suggesting lineage-related molecular functions. Consistent with a key role in human MSC proliferation and survival, we demonstrate that knockdown of MIR181A1HG in the absence of osteogenic stimuli impedes cell cycle progression. Loss of MIR181A1HG enhances differentiation into osteo-chondroprogenitors that produce multiple extracellular matrix proteins. RNA-seq analysis shows that loss of chromatin-bound MIR181A1HG alters expression and BMP2 responsiveness of skeletal gene networks (e.g., SOX5 and DLX5). We propose that MIR181A1HG is a novel epigenetic regulator of early stages of mesenchymal lineage commitment towards osteo-chondroprogenitors. This discovery permits consideration of MIR181A1HG and its associated regulatory pathways as targets for promoting new bone formation in skeletal disorders.

Titanium with nanotopography attenuates the osteoclast-induced disruption of osteoblast differentiation by regulating histone methylation.

The bone remodeling process is crucial for titanium (Ti) osseointegration and involves the crosstalk between osteoclasts and osteoblasts. Considering the high osteogenic potential of Ti with nanotopography (Ti Nano) and that osteoclasts inhibit osteoblast differentiation, we hypothesized that nanotopography attenuate the osteoclast-induced disruption of osteoblast differentiation. Osteoblasts were co-cultured with osteoclasts on Ti Nano and Ti Control and non-co-cultured osteoblasts were used as control. Gene expression analysis using RNAseq showed that osteoclasts downregulated the expression of osteoblast marker genes and upregulated genes related to histone modification and chromatin organization in osteoblasts grown on both Ti surfaces. Osteoclasts also inhibited the mRNA and protein expression of osteoblast markers, and such effect was attenuated by Ti Nano. Also, osteoclasts increased the protein expression of H3K9me2, H3K27me3 and EZH2 in osteoblasts grown on both Ti surfaces. ChIP assay revealed that osteoclasts increased accumulation of H3K27me3 that represses the promoter regions of Runx2 and Alpl in osteoblasts grown on Ti Control, which was reduced by Ti Nano. In conclusion, these data show that despite osteoclast inhibition of osteoblasts grown on both Ti Control and Ti Nano, the nanotopography attenuates the osteoclast-induced disruption of osteoblast differentiation by preventing the increase of H3K27me3 accumulation that represses the promoter regions of some key osteoblast marker genes. These findings highlight the epigenetic mechanisms triggered by nanotopography to protect osteoblasts from the deleterious effects of osteoclasts, which modulate the process of bone remodeling and may benefit the osseointegration of Ti implants.

Chromatin dynamics regulate mesenchymal stem cell lineage specification and differentiation to osteogenesis.

Multipotent mesenchymal stromal cells (MSCs) are critical for regeneration of multiple tissues. Epigenetic mechanisms are fundamental regulators of lineage specification and cell fate, and as such, we addressed the question of which epigenetic modifications characterize the transition of nascent MSCs to a tissue specific MSC-derived phenotype. By profiling the temporal changes of seven histone marks correlated to gene expression during proliferation, early commitment, matrix deposition, and mineralization stages, we identified distinct epigenetic mechanisms that regulate transcriptional programs necessary for tissue-specific phenotype development. Patterns of stage-specific enrichment of histone modifications revealed distinct modes of repression and activation of gene expression that would not be detected using single endpoint analysis. We discovered that at commitment, H3K27me3 is removed from genes that are upregulated and is not acquired on downregulated genes. Additionally, we found that the absence of H3K4me3 modification at promoters defined a subset of osteoblast-specific upregulated genes, indicating that acquisition of acetyl modifications drive activation of these genes. Significantly, loss or gain of H3K36me3 was the primary predictor of dynamic changes in temporal gene expression. Using unsupervised pattern discovery analysis the signature of osteogenic-related histone modifications identified novel functional cis regulatory modules associated with enhancer regions that control tissue-specific genes. Our work provides a cornerstone to understand the epigenetic regulation of transcriptional programs that are important for MSC lineage commitment and lineage, as well as insights to facilitate MSC-based therapeutic interventions.

Oncofetal Epigenetic Bivalency in Breast Cancer Cells: H3K4 and H3K27 Tri-Methylation as a Biomarker for Phenotypic Plasticity.

Alterations in the epigenetic landscape are fundamental drivers of aberrant gene expression that contribute to cancer progression and pathology. Understanding specific modes of epigenetic regulation can be used to identify novel biomarkers or targets for therapeutic intervention to clinically treat solid tumors and leukemias. The bivalent marking of gene promoters by H3K4me3 and H3K27me3 is a primary mechanism to poise genes for expression in pluripotent embryonic stem cells (ESC). In this study we interrogated three well-established mammary cell lines to model epigenetic programming observed among breast cancer subtypes. Evidence is provided for a distinct bivalent signature, activating and repressive histone marks co-residing at the same gene promoter, in the MCF7 (ESR/PGR+) luminal breast cancer cell line. We identified a subset of genes, enriched for developmental pathways that regulate cellular phenotype and signaling, and partially recapitulate the bivalent character observed in ESC. We validated the biological relevance of this "oncofetal epigenetic" signature using data from ESR/PGR+ tumor samples from breast cancer patients. This signature of oncofetal epigenetic control is an informative biomarker and may provide novel therapeutic targets, selective for both recurring and treatment-resistant cancers. J. Cell. Physiol. 231: 2474-2481, 2016. © 2016 Wiley Periodicals, Inc.

MicroRNA-378-mediated suppression of Runx1 alleviates the aggressive phenotype of triple-negative MDA-MB-231 human breast cancer cells.

The Runx1 transcription factor, known for its essential role in normal hematopoiesis, was reported in limited studies to be mutated or associated with human breast tumor tissues. Runx1 increases concomitantly with disease progression in the MMTV-PyMT transgenic mouse model of breast cancer. Compelling questions relate to mechanisms that regulate Runx1 expression in breast cancer. Here, we tested the hypothesis that dysregulation of Runx1-targeting microRNAs (miRNAs) allows for pathologic increase of Runx1 during breast cancer progression. Microarray profiling of the MMTV-PyMT model revealed significant downregulation of numerous miRNAs predicted to target Runx1. One of these, miR-378, was inversely correlated with Runx1 expression during breast cancer progression in mice and in human breast cancer cell lines MCF7 and triple-negative MDA-MB-231 that represent early- and late-stage diseases, respectively. MiR-378 is nearly absent in MDA-MB-231 cells. Luciferase reporter assays revealed that miR-378 binds the Runx1 3' untranslated region (3'UTR) and inhibits Runx1 expression. Functionally, we demonstrated that ectopic expression of miR-378 in MDA-MB-231 cells inhibited Runx1 and suppressed migration and invasion, while inhibition of miR-378 in MCF7 cells increased Runx1 levels and cell migration. Depletion of Runx1 in late-stage breast cancer cells resulted in increased expression of both the miR-378 host gene PPARGC1B and pre-miR-378, suggesting a feedback loop. Taken together, our study identifies a novel and clinically relevant mechanism for regulation of Runx1 in breast cancer that is mediated by a PPARGC1B-miR-378-Runx1 regulatory pathway. Our results highlight the translational potential of miRNA replacement therapy for inhibiting Runx1 in breast cancer.

Histone H3 lysine 4 acetylation and methylation dynamics define breast cancer subtypes.

The onset and progression of breast cancer are linked to genetic and epigenetic changes that alter the normal programming of cells. Epigenetic modifications of DNA and histones contribute to chromatin structure that result in the activation or repression of gene expression. Several epigenetic pathways have been shown to be highly deregulated in cancer cells. Targeting specific histone modifications represents a viable strategy to prevent oncogenic transformation, tumor growth or metastasis. Methylation of histone H3 lysine 4 has been extensively studied and shown to mark genes for expression; however this residue can also be acetylated and the specific function of this alteration is less well known. To define the relative roles of histone H3 methylation (H3K4me3) and acetylation (H3K4ac) in breast cancer, we determined genomic regions enriched for both marks in normal-like (MCF10A), transformed (MCF7) and metastatic (MDA-MB-231) cells using a genome-wide ChIP-Seq approach. Our data revealed a genome-wide gain of H3K4ac associated with both early and late breast cancer cell phenotypes, while gain of H3K4me3 was predominantly associated with late stage cancer cells. Enrichment of H3K4ac was over-represented at promoters of genes associated with cancer-related phenotypic traits, such as estrogen response and epithelial-to-mesenchymal transition pathways. Our findings highlight an important role for H3K4ac in predicting epigenetic changes associated with early stages of transformation. In addition, our data provide a valuable resource for understanding epigenetic signatures that correlate with known breast cancer-associated oncogenic pathways.

Genome-wide co-occupancy of AML1-ETO and N-CoR defines the t(8;21) AML signature in leukemic cells.

BACKGROUND: Many leukemias result from chromosomal rearrangements. The t(8;21) chromosomal translocation produces AML1-ETO, an oncogenic fusion protein that compromises the function of AML1, a transcription factor critical for myeloid cell differentiation. Because of the pressing need for new therapies in the treatment of acute myleoid leukemia, we investigated the genome-wide occupancy of AML1-ETO in leukemic cells to discover novel regulatory mechanisms involving AML-ETO bound genes. RESULTS: We report the co-localization of AML1-ETO with the N-CoR co-repressor to be primarily on genomic regions distal to transcriptional start sites (TSSs). These regions exhibit over-representation of the motif for PU.1, a key hematopoietic regulator and member of the ETS family of transcription factors. A significant discovery of our study is that genes co-occupied by AML1-ETO and N-CoR (e.g., TYROBP and LAPTM5) are associated with the leukemic phenotype, as determined by analyses of gene ontology and by the observation that these genes are predominantly up-regulated upon AML1-ETO depletion. In contrast, the AML1-ETO/p300 gene network is less responsive to AML1-ETO depletion and less associated with the differentiation block characteristic of leukemic cells. Furthermore, a substantial fraction of AML1-ETO/p300 co-localization occurs near TSSs in promoter regions associated with transcriptionally active loci. CONCLUSIONS: Our findings establish a novel and dominant t(8;21) AML leukemia signature characterized by occupancy of AML1-ETO/N-CoR at promoter-distal genomic regions enriched in motifs for myeloid differentiation factors, thus providing mechanistic insight into the leukemic phenotype.

Chromatin modifiers and histone modifications in bone formation, regeneration, and therapeutic intervention for bone-related disease.

Post-translational modifications of chromatin such as DNA methylation and different types of histone acetylation, methylation and phosphorylation are well-appreciated epigenetic mechanisms that confer information to progeny cells during lineage commitment. These distinct epigenetic modifications have defined roles in bone, development, tissue regeneration, cell commitment and differentiation, as well as disease etiologies. In this review, we discuss the role of these chromatin modifications and the enzymes regulating these marks (methyltransferases, demethylases, acetyltransferases, and deacetylases) in progenitor cells, osteoblasts and bone-related cells. In addition, the clinical relevance of deregulated histone modifications and enzymes as well as current and potential therapeutic interventions targeting chromatin modifiers are addressed.

Could lncRNAs be the missing links in control of mesenchymal stem cell differentiation?

Long suspected, recently recognized, and increasingly studied, non protein-coding RNAs (ncRNAs) are emerging as key drivers of biological control and pathology. Since their discovery in 1993, microRNAs (miRNAs) have been the subject of intense research focus and investigations have revealed striking findings, establishing that these molecules can exert a substantial level of biological control in numerous tissues. More recently, long ncRNAs (lncRNAs), the lesser-studied siblings of miRNA, have been suggested to have a similar robust role in developmental and adult tissue regulation. Mesenchymal stem cells (MSCs) are an important source of multipotent cells for normal and therapeutic tissue repair. Much is known about the critical role of miRNAs in biogenesis and differentiation of MSCs however; recent studies have suggested lncRNAs may play an equally important role in the regulation of these cells. Here we highlight the role of lncRNAs in the regulation of mesenchymal stem cell lineages including adipocytes, chondrocytes, myoblasts, and osteoblasts. In addition, the potential for these noncoding RNAs to be used as biomarkers for disease or therapeutic targets is also discussed.

Epigenetic pathways regulating bone homeostasis: potential targeting for intervention of skeletal disorders.

Epigenetic regulation utilizes different mechanisms to convey heritable traits to progeny cells that are independent of DNA sequence, including DNA silencing, post-translational modifications of histone proteins, and the post-transcriptional modulation of RNA transcript levels by non-coding RNAs. Although long non-coding RNAs have recently emerged as important regulators of gene imprinting, their functions during osteogenesis are as yet unexplored. In contrast, microRNAs (miRNAs) are well characterized for their control of osteogenic and osteoclastic pathways; thus, further defining how gene regulatory networks essential for skeleton functions are coordinated and finely tuned through the activities of miRNAs. Roles of miRNAs are constantly expanding as new studies uncover associations with skeletal disorders. The distinct functions of epigenetic regulators and evidence for integrating their activities to control normal bone gene expression and bone disease will be presented. In addition, potential for using "signature miRNAs" to identify, manage, and therapeutically treat osteosarcoma will be discussed in this review.

The bone-specific Runx2-P1 promoter displays conserved three-dimensional chromatin structure with the syntenic Supt3h promoter.

Three-dimensional organization of chromatin is fundamental for transcriptional regulation. Tissue-specific transcriptional programs are orchestrated by transcription factors and epigenetic regulators. The RUNX2 transcription factor is required for differentiation of precursor cells into mature osteoblasts. Although organization and control of the bone-specific Runx2-P1 promoter have been studied extensively, long-range regulation has not been explored. In this study, we investigated higher-order organization of the Runx2-P1 promoter during osteoblast differentiation. Mining the ENCODE database revealed interactions between Runx2-P1 and Supt3h promoters in several non-mesenchymal human cell lines. Supt3h is a ubiquitously expressed gene located within the first intron of Runx2. These two genes show shared synteny across species from humans to sponges. Chromosome conformation capture analysis in the murine pre-osteoblastic MC3T3-E1 cell line revealed increased contact frequency between Runx2-P1 and Supt3h promoters during differentiation. This increase was accompanied by enhanced DNaseI hypersensitivity along with RUNX2 and CTCF binding at the Supt3h promoter. Furthermore, interplasmid-3C and luciferase reporter assays showed that the Supt3h promoter can modulate Runx2-P1 activity via direct association. Taken together, our data demonstrate physical proximity between Runx2-P1 and Supt3h promoters, consistent with their syntenic nature. Importantly, we identify the Supt3h promoter as a potential regulator of the bone-specific Runx2-P1 promoter.

Epigenetic landscape during osteoblastogenesis defines a differentiation-dependent Runx2 promoter region.

Runx2 is a developmentally regulated gene in vertebrates and is essential for bone formation and skeletal homeostasis. The induction of runx2-P1 isoform transcripts is a hallmark of early osteoblastogenesis. Although previous in vitro studies have defined a minimal Runx2-P1 promoter sequence with well-characterized functional elements, several lines of evidence suggest that transcription of the Runx2-P1 isoform relies on elements that extend beyond the previously defined P1 promoter boundaries. In this study, we examined Runx2-P1 transcriptional regulation in a cellular in vivo context during early osteoblastogenesis of MC3T3-E1 cultures and BMSCs induced towards the bone lineage by multi-layered analysis of the Runx2-P1 gene promoter using the following methodologies: 1) sequence homology among several mammalian species, 2) DNaseI hypersensitivity coupled with massively parallel sequencing (DNase-seq), and 3) chromatin immunoprecipitation of activating histone modifications coupled with massively parallel sequencing (ChIP-seq). These epigenetic features have allowed the demarcation of boundaries that redefine the minimal Runx2-P1 promoter to include a 336-bp sequence that mediates responsiveness to osteoblast differentiation. We also find that an additional level of control is contributed by a regulatory region in the 5'-UTR of Runx2-P1.

Disruption of crosstalk between mesenchymal stromal and tumor cells in bone marrow as a therapeutic target to prevent metastatic bone disease.

Skeletal metastasis is a serious complication of many primary cancers. A common feature of tumor cells that metastasize to the bone marrow microenvironment is that they initiate a cascade of events, recruiting and presumably/potentially altering the phenotype of bone marrow mesenchymal stromal cells (MSC) to produce an environment that allows for tumor growth and in some cases, drug-resistant dormancy of latent cancer cells. Consequently the MSC population can contribute to metastatic disease through several distinct mechanisms by differentiating into cancer-associated fibroblasts (CAFs). Understanding the expression and epigenetic changes that occur as normal MSCs become associated with metastatic tumors would reveal possible therapeutic targets for treating skeletal metastasis.

Genomic occupancy of Runx2 with global expression profiling identifies a novel dimension to control of osteoblastogenesis.

BACKGROUND: Osteogenesis is a highly regulated developmental process and continues during the turnover and repair of mature bone. Runx2, the master regulator of osteoblastogenesis, directs a transcriptional program essential for bone formation through genetic and epigenetic mechanisms. While individual Runx2 gene targets have been identified, further insights into the broad spectrum of Runx2 functions required for osteogenesis are needed. RESULTS: By performing genome-wide characterization of Runx2 binding at the three major stages of osteoblast differentiation--proliferation, matrix deposition and mineralization--we identify Runx2-dependent regulatory networks driving bone formation. Using chromatin immunoprecipitation followed by high-throughput sequencing over the course of these stages, we identify approximately 80,000 significantly enriched regions of Runx2 binding throughout the mouse genome. These binding events exhibit distinct patterns during osteogenesis, and are associated with proximal promoters and also non-promoter regions: upstream, introns, exons, transcription termination site regions, and intergenic regions. These peaks were partitioned into clusters that are associated with genes in complex biological processes that support bone formation. Using Affymetrix expression profiling of differentiating osteoblasts depleted of Runx2, we identify novel Runx2 targets including Ezh2, a critical epigenetic regulator; Crabp2, a retinoic acid signaling component; Adamts4 and Tnfrsf19, two remodelers of the extracellular matrix. We demonstrate by luciferase assays that these novel biological targets are regulated by Runx2 occupancy at non-promoter regions. CONCLUSIONS: Our data establish that Runx2 interactions with chromatin across the genome reveal novel genes, pathways and transcriptional mechanisms that contribute to the regulation of osteoblastogenesis.